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Abstract 

A technique for extracting from the appropriate field equations the 
relativistic motion of Schwarzschild, Reissner-Nordstrom and Kerr par¬ 
ticles moving in external fields is motivated and illustrated. The key 
assumptions are that (a) the particles are isolated and (b) near the 
particles the wave fronts of the radiation generated by their motion 
are smoothly deformed spheres. No divergent integrals arise in this 
approach. The particles are not test particles. The formalism is used, 
however, to derive the Mathisson-Papapetrou equations of motion of 
spinning test particles, neglecting spin-spin terms. 


1 Introduction 

We describe a method for modeling the relativistic motion of (uncharged 
or charged) particles moving in external gravitational and electromagnetic 
fields in general relativity. The particles modify the fields in which they 
move and thus are not test particles. The modified fields near the particles 
are predominantly the Schwarschild, Reissner-Nordstrom or Kerr fields and 
thus the particles will be referred to as Schwarschild, Reissner-Nordstrom 
or Kerr particles respectively. The origin of the approach described here is 
the seminal and challenging paper by Robinson and Robinson [I], with the 
disconcertingly modest title of “Equations of Motion in the Linear Approxi¬ 
mation”, published in the festschrift in honor of Professor J. L. Synge. This 
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paper inspired early work involving space-times which are more special than 
those required here and which can be found in a series of papers by Hogan 
and Imaeda (see also 0). A line element, which is a byproduct of 

a study of gravitational radiation from bounded sources [6], plays a key role 
in the development of our method. An early application can be found in [7]. 
The most up-to-date references to our technique and applications are in the 
paper [8] and the books [9] and m- Further development of this work has 
also involved graduate students Takashi Fukumoto and Shinpei Ogawa in 
Tohoku University, Aishling Nic an Tuile in University College Dublin and 
Florian Bolgar in Ecole Normale Superieure, Paris and has taken the form 
of M.Sc. and Ph.D. theses and an Internship Report. 

To introduce our approach we begin with the Eddington-Finkelstein 
form of the Schwarzschild line element: 

+ '^2/^) + 2 dn dr -b -^ du^ , (1) 

with 

Po = 1 + + 2/^) ■ (2) 

Here x, y are stereographic coordinates taking values in the ranges — oo < 
X < -boo and —oo<y< -boo, u is a null coordinate (in the sense that the 
hypersurfaces u = constant are null) with —oo < u < -boo. The generators 
of these null hypersurfaces are null geodesics labelled by {x, y) and r (with 
0 < r < -boo) is an affine parameter along them. When the mass parameter 
m = 0 the space-time is Minkowskian with line element 

dsg = —r^pQ ^(dx^ -b dy^) + 2dudr + du^ = ytjdX^dX^ , (3) 


with A* = (A, Y, Z, T) rectangular Cartesian coordinates and time and ryj = 
diag(—1, —1, —1,-bl). In this background space-time with line element ([3|) 
the hypersurfaces u = constant are future null cones with vertices on the 
time-like geodesic r = 0. We can consider 


, i , n 2 m 2 

jij dx dx-' = - du 


(4) 


in coordinates x* = {x,y,r,u), as a perturbation of this background which 
is singular on the time-like geodesic r = 0 and which, when added to the 
line element ([3|) produces the line element ([T|). Although we have no need 
to assume it now, we shall, in the sequel, assume that m is small of first 
order and write m = Oi in order to simplify the complicated calculations 
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that develop. We shall consider the background space-time a model of 
the external field in which the small mass is located. In the example here 
the background space-time is Minkowskian space-time and thus there is no 
external field present. Also in this example the time-like geodesic equations 
(satisfied by r = 0 in ([3])) are the equations of motion of the small mass. 

In the space-time with line element (l3|) let A* = w'^{u) be an arbitrary 
time-like world line with u proper time or arc length along it. Then r>*(n) = 
dw^/du is the unit tangent to this line and satisfies VjV^ = rjijv^v^ = +1 
along the world line. Equivalently n* is the 4-velocity of the particle with 
world line A* = w'^{u). Thus a* = dv^/du is the 4-acceleration and satisfies 
OjV^ = r]ijadv^ = 0 along the world line. We shall lower and raise indices, 
referring to the coordinates A*, using r]ij and respectively, with the latter 
defined by rjjk = ■ We note that n® and ad are the components of the 4- 

velocity and 4-acceleration of the particle expressed in the coordinates A®. 
These features of the world line can be incorporated in the Minkowskian 
line element by replacing ([3|) by 

dso = —+ dy^) + 2dudr + {1 — 2hor)du^ , (5) 


with 

Po = xn^(u)+ ?/n^(n) + |l - + y^)| , 

( 6 ) 

and 

ho = ^(log Po) = Oj , (7) 

with 


-x5{ - 7/5^ - |l - i(x^ + y^)| (5^ + |l + ^(x^ + y^)| 5l . 

( 8 ) 

In ([5]) r = 0 is the world line A® = w^{u) with the 4-velocity components 
appearing in ([5]) via the function Pq and the 4-acceleration appearing via /iq. 
The future null cones with vertices on this world line are the null hypersur¬ 
faces u = constant. The generators of these null cones are the null geodesic 
integral curves of the vector field , which satisfies kjk^ = 0 and kjV^ = -|-1. 
On each future null cone the generators are labeled by the stereographic co¬ 
ordinates X, y and r > 0 is an affine parameter along them. To see how 
Einstein’s vacuum held equations can lead to the equations of motion as we 
have dehned them in the paragraph above we consider a space-time with 


y = Pq-i 


3 




line element 


= —t^Pq ^(dx^ + dy^) + 2dudr + — 2hor — ^ (^g) 

which is a generalization of the Schwarzschild space-time with line element 
©• The Ricci tensor components in coordinates x* = {x,y,r,u) are now 
given exactly via 

Rijdx^dx^ = _ (^ 10 ) 

With m 7 ^ 0 we see that the vacuum field equations Rij = 0 require ajk^ = 
ho = 0 and this should hold for all possible kR Hence we must have, as 
equations of motion, the time-like geodesic equations 

= 0 , ( 11 ) 

for the world line r = 0 in the Minkowskian background space-time. This 
simple example illustrates a connection between the equations of motion 
and the field equations. In general we shall find, however, that the field 
equations by themselves will not be sufficient to determine the equations 
of motion. They will have to be supplemented by a smoothness condition 
on the wave fronts of the possible radiation produced by the motion of the 
particle. A singular point on a wave front translates in space-time into a 
singularity along a generator, extending to future null infinity, of the null 
hypersurface history of the wave front and thus violating the concept of an 
isolated source. A useful example with which to illustrate this is to consider 
the relativistic motion of a Reissner-Nordstrom particle moving in external 
electromagnetic and gravitational fields. 

2 Reissner-Nordstrom Particle 

We construct here a model of a Reissner-Nordstrom particle of small mass 
m = Oi and small charge e = Oi moving in external vacuum gravita¬ 
tional and electromagnetic fields. In this case the background space-time 
is a general solution of the vacuum Einstein-Maxwell field equations. We 
shall require a knowledge of this space-time, and the Maxwell field, in the 
neighborhood of an arbitrary time-like world line r = 0 (say). The charged 
mass will then be introduced as a perturbation of this space-time, which is 
singular on this world line in the background space-time (in similar fash¬ 
ion to dH) in the space-time with line element ([3|)), in such a way that the 
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perturbed space-time is, for small values of r, predominantly the Reissner- 
Nordstrom space-time and Maxwell field. Since we have indicated in Section 
1 that we shall make use of the field equations along with conditions on the 
wave fronts of possible radiation (electromagnetic and/or gravitational in 
the current example) produced by the motion of our particle, and since the 
histories of wave fronts are null hypersurfaces in space-time, we shall work 
in the background space-time and in the perturbed background space-time 
in a coordinate system based on a family of null hypersurfaces. Thus we 
begin by writing the line element of the background space-time as 

ds^ = - (^^2)2 = gab^'^ , (12) 

with the 1-forms and given by 

= rp“^(e" cosh/3 dx + sinh/3 dy + adu) , (13) 

7^2 _ rp“^(e" sinh/3dx + cosh/3 dy + 6du) , (14) 

= dr + ^ dn , (15) 

= du . (16) 

The constants gab = dba constitute the components of the metric tensor 
on the tetrad defined by the basis 1-forms. Tetrad indices will be denoted 
by the early letters of the alphabet a, 6, c,... , and tetrad indices will be 
lowered with gab and raised with where g°'^gbc = d/- This line element is 
completely general in that it involves six functions p, a, /3, a, b, c each of the 
four coordinates x, y, r, u. It is thus equivalent to line elements constructed 
by Sachs m and by Newman and Unti m- It was constructed in [6] for the 
study of gravitational radiation from bounded sources and was specifically 
designed to have the Robinson-Trautman m, M) line elements appear 
as a convenient special case (got by simply putting a = f3 = 0). This latter 
property also makes it convenient in the current context. We shall take 
r = 0 to be the equation of an arbitrary time-like world line in the space- 
time with line element (m. The space-time in the neighborhood of this 
world line has the well-known Fermi property (see, for example, m, US]) 
that there exist coordinates in terms of which the metric tensor components 
satisfy 

9ij = Vij + O(r^) . (17) 

In view of the line element ([5]) we can implement this property in the context 
of the line element (jl2j) by expanding the six functions appearing there in 
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powers of r as follows: 

P = Po{l + q2r‘^ + + ■ ■ ■ ), (18) 

a = a2 + as + ... , (19) 

/3 = + ... , (20) 

a = air + a2r^ + ... , (21) 

b = hir + h2r‘^ + ... , (22) 

c = Co + Cl r + C2 + ... , (23) 


where the coefficients of the powers of r are functions of x, y, u. Following 
(l5|) the function Pq here is given by ([6]) and 

co = l = AlogPo with A = P^(^-^ + ^^ , (24) 

and 

Cl = -2ho , (25) 

with ho given by ([7|). The operator A is the Laplacian on the unit 2-sphere. 
The relationship between the rectangular Cartesian coordinates and time A* 
and the coordinates x,y,r,u near the arbitrary time-like world line r = 0 
(<t4> A* = w'^{u), introduced in Section 1) is given by (see, for example, [T7] i 

= w\u) + rk\ (26) 


neglecting 0(r^)-terms. 

At this point in the discussion we should note some formulas associated 
with Minkowskian space-time with line element ([5]) or, equivalently, with 
metric tensor components r]ij in coordinates A*. A useful basis of 1-forms 
suggested by ([5]) is 


= r Pq ^dx , u? = r Pq ^dy , = dr + ^ — ho du , = du . (27) 

When expressed in terms of the coordinates A* using the transformation 
(|^ these read 


= -Po ^ dX^ = -oji , 
ox 

dy 

= (^i - dX'- = W4 , 

= ki dA* = ws . 


(28) 

(29) 

(30) 

(31) 
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In coordinates X* the vector fields with components ,dk'^/dx^dk^ /dy 

form a basis and it is useful to express the metric tensor components on 
this basis as 


rii — _ p2 
7 — 


dy dy dy 

dx dx dy dy 


+ yy + yy-yy 


(32) 


Also the second partial derivatives of y play a key role in the calculations 
and when they are expressed on this basis we have 


d'^y 

dx"^ 


(33) 

d'^y 


(34) 

d'^y 

dxdy 


(35) 

As a simple indication of the usefulness of these formulas we notice 
we add pSP and (l3H) we obtain 

that if 


Ay + 2y = 2y , 

(36) 

with the operator A defined in (I24h. Taking the scalar product, with re¬ 
spect to the Minkowskian metric tensor r]ij, of this equation with the 4- 

acceleration ad 

yields 



A/iq -|- 2 ho = 0 , 

(37) 


with ho given by This demonstrates that /iq is an f = 1 spherical 

harmonic. A spherical harmonic Q of order I is a smooth solution, for 
—oo < x,y < + 00 , of the equation AQ + l{l + 1)Q = 0 for f = 0,1,2,... . 

Returning now to the construction of the background Einstein-Maxwell 
field, the potential l~form leading to the background electromagnetic field 
is given by 

A = L dx + M dy + K du , (38) 

where L, M, K are functions of x,y,r,u. We have used the freedom to add 
an exact differential (i.e. a gauge transformation) to remove one component 
of this 1-form for convenience. The functions appearing in (j38p will be 
expanded in powers of r in such a way that the resulting Maxwell 2-form 
F = dA, the exterior derivative of (I38p . is non-singular at r = 0, since this 
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will be the external Maxwell field. The simplest expansions to achieve this 
are: 


L 

= L 2 + T 3 + • • • 

(39) 

M 

= M 2 + M 3 + ... , 

(40) 

K 

= r Xi + X 2 + ... , 

(41) 


with the coefficients of the powers of r functions of x, y, u. The Maxwell 
2 -form is given by 

T = = (42) 

where F^b = —Fba are the components of the 2 -form on the tetrad defined 
by the basis 1-forms (|13p - ()16p . The leading terms, in powers of r, of the 
tetrad components Fab which will interest us are, in light of the expansions 
([IHP-dlSl) and (13^-(1TT]1. 

Fi 3 =-2Po T 2 + 0(r) , F 23 = -2Po M 2 + 0(r) , F 34 = iiTi + 0(r) . (43) 


Thus on r = 0 we can replace the basis d-CSI) by the Minkowskian 
counterparts ([28]) - (l3Tp and the coordinates x, y, r, u by the coordinates X* 
following (f26j) to arrive at 


1 • 31^^ 1 • 

L, = -F,,{u)F—,M, = -F,,{u)F — 


K,=Fi,{u)Fv^ , (44) 


where Fij{u) = —Fji{u) are the components of the external Maxwell tensor 
on the world line r = 0 calculated in the coordinates X*. With F given 
by dS]) above we see that now the dependence of the functions L 2 ,M 2 ,Ki 
on X, y is explicitly known. We can therefore use the useful formulas ([33])- 
(I35p to verify that Maxwell’s vacuum held equations, d*F = 0, where the 
star denotes the Hodge dual of the 2-form F, are satished by L 2 , M 2 , Ki. 
The relevant equations obtained from the leading terms, in powers of r, in 
Maxwell’s equations are 


Xi = 


dL2 ^ dM2 


and 


dy 
+ 2 L 2 


, AXi + 2H[ 


dx 

dKi 
dx 


dL2 ^ dM2 


dy \ 
d 




dx dy 

dM2 dL2 


= 0 , 


dx dy 

dM2 dL2 


dx 


dy 


= 0 , 
= 0 


(45) 

(46) 

(47) 
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We see that the second equation in (HSl) is a consequence of ()^6]) and (|^71) 
and the two equations in (j45p imply that Ki is an / = 1 spherical harmonic 
since 

AKi+2Ki = 0 . (48) 

This latter equation can, of course, be verified directly using the last of (1141) 
and the formulas (f33]l - (f35]) . 

The analogue, for the gravitational field, of the Maxwell 2-form with 
tetrad components Fab given by ()43p is the Weyl conformal curvature tensor 
with tetrad components Cabcd- The components which will be of interest to 
us in the neighborhood of the world line r = 0 are given by 

^1313 + *(^1323 = 6 ( 0:2 + f/ 32 ) + 0{r) , (49) 

and 

C 3431 + fC3432 = |ai + ibi + 4p2 + 0(r) , (50) 

where ( = x + iy and a bar denotes complex conjugation. These must be 
examined in conjunction with the Einstein-Maxwell vacuum field equations 
for the background space-time: 

Rab = 2 i^FcaF^^b - -^gabFcd.F'^^^ = 2 Eab , (51) 

where Rab are the tetrad components of the Ricci tensor of the background 
space-time, gab are the tetrad components of the metric tensor and, as 
indicated following (1161) . tetrad indices are raised using the components 
of the inverse of the matrix with components gab and Eab = Fba are the 
tetrad components of the electromagnetic energy-momentum tensor. To 
satisfy the field equation R 33 = 2 E 33 + 0{r) we find that 

q2 = ^P^{Ll+Mi) . (52) 

This can be rewritten, using L 2 and M 2 given by (j44|) and the expression 
(I32p for the components of the inverse of the Minkowskian metric tensor, 
simply as 

q2 = Fpj{u) P P . (53) 

Here Fpj[u) are the components of the Maxwell tensor on the world line 
r = 0 in coordinates X* while FPi(u) = gPMji[u). Using this expression for 
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q 2 in (l50l) and the basis ([28]) - (l3T]) . in similar fashion to the passage from 
(li3]l to ()i^ above, we conclude from and ([501) that 

1 Bhi f)h^ 1 Bhi Bh^ 

a2 = -^Pi ^2 = -^Pi (54) 

and 

Oi = -Pq \ Cijki{u) P k’" — + F^i{u) Fpj{u) P —j , (55) 

hi = -Pi [Cijui{u)PPk'^^ + F^^{u)Fp,{u)P—y (56) 

where Cijki{u) are the components of the Weyl tensor of the background 
space-time calculated on the world line r = 0 in the coordinates X*. With 
Pq given by ([6]) and by ([8]) , we see that the dependence of q 2 in (1531) , a 2 and 
h in (fMl) . and oi and hi in (l5^ and (1561) on the coordinates x, y is explicitly 
known. The Einstein-Maxwell vacuum field equations R 12 = 2 Ei 2 + 0{r) 
and Rii — R 22 = — E 22 ) + 0 {r) yield the pair of real field equations 

incorporated in the single complex equation: 

2{a2 + i/ 32 ) =^cbi + ibi + 4 pi■ (57) 

The pair of real field equations Ris = 2 E 13 + 0{r) and R 23 = 2 E 23 + 0{r) 
can be written as the complex equation; 

ai + ib, + 4P2 % = 2 Pi ^{Po^{a 2 + i/Ja)} • (58) 

The field equations (1571) and (l58]l must be satisfied by (72 j «2 , /32 , oi, bi given 
by (f53]l " (f56]l . This important check can be carried out using the useful 
formulas (p3]l - (f35]) . We are now ready to introduce the Reissner-Nordstrom 
particle as a perturbation of this background space-time and Maxwell field. 


The line element of the perturbed space-time is given by 

ds^ = -(^91)2 - , (59) 

with the 1-forms , 13“^, and r given by 

^9^ = rp“^(e" cosh/3 dx + e“" sinh/3 dy + adn) , (60) 

1)2 = rp“^(e" sinh/3dx + e““ cosh/3 dy + 6du) , (61) 

= dr + ^du , (62) 

= du . (63) 
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To achieve our aim of having a predominantly Reissner-Nordstrom field for 
small values of r we take the functions appearing here to have the following 
expansions in powers of r: 


p = 

Po{l + hk + q:ik + ... ) , 

(64) 

a = 

012 ^ + Oi^k + ... , 

(65) 

^ = 

hk + + ■ ■ ■ , 

( 66 ) 

a = 

^ — 1 ^2 

-h ao + ai r + a 2 r + ... , 

r 

(67) 

h = 

+ 60 + ^ + &2 ^^ + • • • , 

r 

( 68 ) 

c = 

(P 2 (m-b 2 /_i) , . , . , . 2 , 

^2 j. -t- Co -b Cl r -t- C 2 r -b ... 

(69) 


The coefficients of the various powers of r here are functions of x,y,u. The 
hatted functions differ from their background values by Oi-terms. Thus 


in particular a_i = Oi,ao = Oi,6_i = Oi,6o = Oi, but /_i = O 2 . The 
perturbed potential 1-form 

A = L dx + M dy + K du , (70) 

is predominantly the Lienard-Wiechert 1-form (= e (r“^ — ho) du) up to a 
gauge term so that 

L = r‘^L 2 + r'^U + ... , (71) 

M = Tr^Ms + ... (72) 

k = (g + -^-i) _e/to + rRi+r^i^2 + ... , (73) 

r 


where again the coefficients of the powers of r here are functions of x,y,u 
which differ from their background values by Oi-terms and -R_i = Oa- The 
expansions (|6^ - (|6^ and (r7T]l - (|73]l appear to us to be the minimal assump¬ 
tions necessary to have a small charged mass with gravitational and elec¬ 
tromagnetic fields predominantly those of a Reissner-Nordstrom particle for 
small values of r. They certainly restrict the model of such a particle moving 
in external gravitational and electromagnetic fields and their generalization 
or otherwise is a topic for further study. 

When the metric tensor given by (l5^ - (l63|) and the potential 1-form (f70]l . 
together with the expansions (IM]l - (l69|) and d7T|) - d73]) . are substituted into 
Maxwell’s vacuum field equations and the Einstein-Maxwell vacuum field 
equations, and all terms are gathered on the left hand side of each equation 
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with zero on the right hand side, we find a finite number of terms involving 
inverse powers of r and an infinite number of terms involving positive (or 
zero) powers of r. In an ideal world we would equate the coefficients of each 
of these powers to zero and from the resulting equations derive the coeffi¬ 
cients in the expansions and (l7T]l - d75|) . However the reality is that 

we can at best make the coefficients small in terms of m and e. How small 
depends upon how accurately we require the equations of motion. Also the 
number of coefficients we require in the expansions (|64p - (|69p and (I7ip ^ (|73l) 
depends upon how accurately we require the equations of motion. For the 
purpose of the present illustration we shall calculate the equations of motion 
in first approximation and thus with an 02“error. The null hypersurfaces 
u = constant in the space-time with line element (|59l) are the histories of 
possible wave fronts of radiation (electromagnetic and/or gravitational) pro¬ 
duced by the motion of the Reissner-Nordstrom particle. For small values of 
r the degenerate metric on these hypersurfaces is given by the line element 

dsl = -r^pQ‘^{dx‘^ + d-ip) , (74) 

where Pq differs from its background value Pq given in ([6l) by Oi-terms and 
thus can be written 

Po = Po{l + Qi+ O 2 ) , (75) 

with Qi(x,y,u) = Oi. We will assume that the line elements (|74p are 
smooth deformations of the line element of a 2-sphere. This means that Qi 
is a well behaved function for —00 < x,y < -|-oo and thus free of singularities 
in X, y. Now the field equations yield 


a_i = 

-4eP^L2 + 02 = 02 , 

(76) 

b-i = 

-4eP^M2 + 02 = 0 i , 

(77) 

Co = 

1 -|- AQi 2 Qi -|- 8 e Fij P O 2 , 

(78) 

^A(AQi 

+ 2 Qi) = 6 m ai P — G e Fij P + O 2 ■ 

(79) 


Both terms on the right hand side of (|79l) are I = 1 spherical harmonics 
and thus (|7^ can be easily integrated (discarding the solution of the I = 0 
spherical harmonic equation which is singular for infinite values of x and y) 
to read 

AQi + 2 Qi = Gma-iP — 6 e Fij P + A{u) + O 2 , (80) 

where A{u) = Oi is arbitrary (an I = 0 spherical harmonic). We note that 
spherical harmonics corresponding to / = 0 or / = 1 in the perturbation Qi 
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in ()75p are trivial, in the sense that the 2-sphere remains a 2-sphere under 
these perturbations, and can be neglected. Since the first two terms in (fSOl) 
are both I = 1 spherical harmonics the solution Qi will necessarily have a 
singularity in x, y (a directional singularity) unless they combine to vanish 
or be at most small of second order. Thus we must have 

maik'’ = e Fij F + O 2 , (81) 

for all possible values of F and thus 

mai = e Fij + O 2 ■ (82) 

We have arrived at the Lorentz equations of motion in first approximation. 

The example above serves to illustrate our approach. If it is continued 
to the next order of approximation the equations of motion which emerge 
are [ 8 ] 

mai = e Fij F^u Fpj hf d’^ + + O 3 , (83) 

where h) = 5) — u* Vj is the projection tensor and the dot denotes differenti¬ 
ation with respect to proper time u. The second term on the right hand side 
here is an 02 -correction to the external 4-force. The third term on the right 
hand side is the electromagnetic radiation reaction 4-force. The fourth term 
on the right hand side is a tail term given by an integral with respect to 
proper time u from —00 to the current proper time u of a vector whose com¬ 
ponents are proportional to the external Maxwell field Fij(u) and involve 
functions of u of integration and a pair of space-like vectors, which occur 
naturally, defined along the world line r = 0 in the background space-time. 
When e = 0 we obtain the equations of motion of a Schwarzschild parti¬ 
cle in second approximation given by the approximate time-like geodesic 
equations 

mai = Os . (84) 

There is a formal similarity between the equations of motion (|83p and equa¬ 
tions of motion derived by DeWitt and Brehme |18] . However DeWitt and 
Brehme have removed an infinite term from their equations of motion and 
have considered a scenario (the world line of a charged test particle in a 
curved space-time, and hence utilize only Maxwell’s equations on a curved 
space-time) which is different to that constructed here. 

Finally, so that the reader is left in no doubt as to the procedure involved 
in the passage from (f80l) to ([82]) . we can give a simplified example of the 
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procedure. If A = cos 6* for 0 < 0 < vr (or equivalently for — 1 < A < +1) 
consider a function / = /(A) satisfying the inhomogeneous I = 1 Legendre 
equation (i.e. the inhomogeneous I = 1 spherical harmonic equation satished 
by a function of the single variable A) with an / = 1 Legendre polynomial 
on the right hand side: 

A/ + 2/ = J-|(l-A2)^|+2/ = tiPi(A) , (85) 

where ti is a constant. For this equation to have a solution which is non¬ 
singular at A = ±1 (which are directional singularities) we must have ti = 0 
(analogous to (18111 ). This easily follows from consideration of the general 
solution 

fW = “F + ^2 A(A) + ts Qi(A) , (86) 

0 

where t 2 ,t 3 are constants and (5i(A) is the Legendre function corresponding 
to / = 1. Of course the required directional singularity-free solution also 
requires us to take ts = 0 here. 

3 Kerr Particle 

We outline here the application of our technique to the construction of a 
model of a spinning particle, or Kerr particle, moving in an external vacuum 
gravitational field, and the derivation of its equations of motion. This work, 
carried out by Shinpei Ogawa and the author, is reported in [TO]. Although 
we consider only the equations of motion in first approximation here, the 
calculations are more intricate than those necessary for the study of the 
Reissner-Nordstrom particle. Consequently we must refer the interested 
reader to m for further details. 

The fundamental aspects of our technique are the construction of a back¬ 
ground space-time with an arbitrary time-like world line in the neighbor¬ 
hood of which the space-time is flat, following from the Fermi property, and 
then the introduction of the particle of interest as a perturbation of this 
space-time. In the flat neighborhood of the world line in the background 
space-time we wrote the line element in the form ([5]) which introduces the 
4-velocity and 4-acceleration of the world line into the line element. To 
consider a particle with spin in this context we introduce the spin variables 
(what will become the three independent components of the angular mo¬ 
mentum per unit mass of the particle) into the line element of Minkowskian 
space-time to accompany the 4-velocity and 4-acceleration already present. 
We do this in such a way that if a* = 0 then the form of the line-element 
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of Minkowskian space-time coincides with the Kerr line element, with three 
components of angular momentum per unit mass, in the special case in 
which the mass m = 0. The latter form of the Kerr line element can be 
found in ([H], p.37). Let s*(n) be the components of the spin vector in 
coordinates X*, defined along the arbitrary time-like world line = w'‘{u) 
to be everywhere orthogonal to the 4-velocity u*(u) and Fermi transported, 
thus. 


s* Uj = 0 and 


ds^ 


= — Oi S-' s 


du '' ^ 

respectively. The equivalent spin tensor is defined by 


(87) 


Sij — ^ijkl 


V S = —s 


Jl 5 


( 88 ) 


where Cijki are the components of the Levi-Civita permutation tensor in 
coordinates X*. We introduce s* into the Minkowskian line element by 
replacing ([261) by 

’ ( 89 ) 

where 

F = s^ ki , (90) 

and the subscripts on F denote partial differentiation with respect to x and 
y. We note that Pq and F are given by ([6]) and ([8]). By (j87l) and the useful 

formulas ()33p - (j35p we easily find that K is an Z = 1 spherical harmonic. 

Also we can rewrite (j89p in the general form of (I26p as 

A* = w^{u) + r K* with K* = H— kj , (91) 


showing that for large r in Minkowskian space-time AT* only differs from 
F by an infinitesimal Lorentz transformation generated by the spin tensor. 
The Minkowskian line element is given by dsg = rjij dX^ dX^ with A* given 
now by (I89p . When A* = rc*(ri) (no longer corresponding to r = 0) is a 
geodesic we can choose u* = b\ and the resulting line element is 


dso = —(r^ + F^)pq “^{dx^ + dy'^) + 2dS |dr + ^(dn 
with po given by ([2]) and 



dX = du — Fy dx + Fxdy . 


(93) 
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We note that the latter is not an exact differential unless the spin van¬ 
ishes. The Kerr solution with the three components of angular momentum 
{m ,m ,m s^) (since now = 0 ) is 

■ ( 94 ) 

If the world line = re* (u) is not a time-like geodesic then the Minkowskian 
line element is algebraically more complicated than (1921) . For simplicity we 
shall henceforth neglect spin-spin terms and assume that the 4-acceleration 
is proportional to the spin on the basis that we expect geodesic motion if the 
spin vanishes. Now the Minkowskian part of the background line element, 
when the world-line = w^{u) is not necessarily a time-like geodesic, is 
given by 

9 9 - 9(9 9 2Pn F,, 2Pn Ft 1 

dsl = —F Pq < dx^ + dy^ H- dx du - 2 — dydu \ +2 dr dX 

+2 Fy dudx — 2 Ft dudy + {I — 2ho r) du^ . (95) 

With this preparation we write the line element of the general background 
space-time as 

ds^ = —(e" cosh/3(ix + sinh/3(iy-k 

+(e" sinh (3 dx + e“" cosh/3 dy + bdXF)^ \ +2dr dX + cdX^ , 




(96) 

p = 

Po (1 + 92 + 93 ■ ) , 

(97) 

a = 

air + a2r‘^ + ... , 

(98) 

/3 = 

/3i r + /32 + ... 

(99) 

a = 

Pq Fy 0-1 

2 + -\- aQ ai T . , 

(100) 

h = 

PqFx b_i , 

2 + -|- 5o T ^ T • • • ) 

(101) 

c = 

Co + Cl r + ... 

(102) 


We now require this background line element to be a solution of Einstein’s 
vacuum field equations 

Rab = ^ ■ (103) 
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Equating to zero the powers of r in the Ricci tensor components Rab we 
find the following expressions for the coefficients of the powers of r in the 
expansions (IWll - (ll02l) : 


qi = 

0 = g2 , ai = 4F/32 , Pi = —4Fa2 , 

(104) 

Ot2 

1 rlh^ 

(105) 

h 


(106) 

Also a_i = 0 = b-i 

and 


ao = —7Poia2Fy — /32Fx) , bo = —7Po ( 02 X 3 : + P2Fy) , 

(107) 

ai 

2 BF 

= 3^0 RijM + 0{F) , 

(108) 

bi 

2 BF 

= jPo RiikM k'v’k‘‘ —+ 0(F) , 

(109) 

and Co = 1 while 



Cl = -2ho- FP^ 

+ 5(oir, - hF.). 

( 110 ) 


The 0(-F)-terms not calculated here will get multiplied by F, as in (|104jl . 
(11071) and (jllOp and will therefore not contribute to the end result since we 
are systematically neglecting spin-spin terms. In these formulas Rijki{u) 
are the components of the Riemann curvature tensor in coordinates X* cal¬ 
culated on the world-line X* = w^{u). As well as satisfying the algebraic 
symmetries of the Riemann tensor they satisfy Einstein’s vacuum field equa¬ 
tions rf^Rijki{u) = 0 on the world -line X* = w^{u). 

We introduce the Kerr particle of small mass m = Oi as a perturbation, 
which is predominantly the Kerr field for small values of r, of this background 
vacuum space-time. The simplest way to achieve this appears to be the 
expansions (remembering that we make no claim to uniqueness; see the 
comments following (1731) above): 

ds'^ = —(e" cosh/3dx-|-e““ sinh/3(i?/-b 

-|-(e“ sinh ^ dx + e““ cosh ^ dy + h dS)^ | + 2dr dTi + c dS^ , 

( 111 ) 
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with 


Po{l + q2r‘^ Fqsr^ + ■■■ ) , 

( 112 ) 

di r + d 2 + ... , 

(113) 

^ir +P 2 P+ ■■■ , 

(114) 

Pq Fy 0_1 

9 + + Oq + Cll c + . . . , 

(115) 

PoFx b_i j ~ 

9 + + 60 T ^ T • • • ) 

(116) 

2 m 

-h Co + Cl r + ... 

r 

(117) 


As always the hatted functions of x, y, u differ from their background values 
by Oi-terms. In addition Pq again has the form (I75p involving the first 
order function Qi{x, y, u). Now we must impose the vacuum field equations 
to be satisfied by the metric given via the line element (jllll) with sufficient 
accuracy to enable us to derive the equations of motion for the world line 
A* = w'‘{u) in the background space-time, in first approximation. The 
vacuum field equations provide us with a set of equations which parallel (17^ - 
(l8T]l in the Reissner-Nordstrom example. In addition to qi = 02,q2 = O 2 
we find that 


do = -7PQ(a2Py -/32 Px) - 2mai + O(mF) + O 2 , (118) 

bo = - 7 PQ(a 2 Px + / 32 Fy) - 2mbi + 0{mF) + O 2 , (119) 

a_i = —AmF bi + O 2 , 6-1 = AmF ai + O 2 , (120) 

and 

Co = 1 + AQi + 2Qi- lOmFPo 

- 8 m {aiFy - biFx;) + O 2 ■ ( 121 ) 

In the ramaining field equation, which yields the analogue of (1791) and thence 
the analogue of (fSnjl . the functions oq and bo are multiplied by m and so 
the uncalculated terms 0{mF) in (I118p and (I119p will not contribute to the 
equations of motion in first approximation (i.e. neglecting 02 “terms). The 
equation we obtain, which is the analogue of (IHUp . is 

(A + 2){Qi + 2m Ji + m J 2 ) = 6m Ui F — 6m Rijki{u) F + A{u) + O 2 , 

( 122 ) 
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where A{u) = Oi is a function of integration as in (IHOl) and Ji and J 2 are 
spin-curvature terms given by 


Ji = P^F R,,ki{u) P v^-^— + -R^Mu) P , 


which is an Z = 3 spherical harmonic and 


J 2 


Po Rijkl{u) Sm t’* 


( dk^ dk'^ 
I dx dy 


dk^ dk'^ 
dy dx 


(123) 


(124) 


which is an Z = 2 spherical harmonic, results which are obtained using the 
useful formulas (j33p - ()35p . The hrst two terms on the right hand side of 
(11221) are I = 1 spherical harmonics and so, just as in the case of (IHop . the 
solution Qi will be free of directional singularities provided the sum of these 
terms is zero or at most small of second order, i.e. 

muiP = mRijki{u) P + O 2 , (125) 

for all possible P and thus we arrive at the equations of motion of the Kerr 
particle in first approximation: 


mai = m Rijki{u) + O 2 ■ (126) 

The algebraic form of the right hand side of these equations is, perhaps, 
what one would expect. The numerical factor of unity distinguishes these 
equations of motion from those of Mathisson-Papapetrou for a spinning test 
particle. We have here a model of a Kerr particle and it is manifestly not 
a test particle. The rest mass m is already small and making it smaller 
does not make it a test particle since so long as m 7^ 0 its presence perturbs 
the background space-time (|96ll . However in the present context we can say 
something about spinning test particles and this is the topic we turn to now. 

4 Spinning Test Particle 

Since we are now concerned with test particles moving in vacuum gravita¬ 
tional fields we will concentrate on what we called background space-times 
above. We hrst consider the non-spinning test particle whose world line, ac¬ 
cording the geodesic hypothesis, is a time-like geodesic. With the geodesic 
hypothesis accepted for a non-spinning test particle it is relatively easy to 
devise a strategy, using the formalism of this paper, which will lead us to the 
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equations of motion of a spinning test particle. As in Section 3 we will neglect 
spin-spin terms, but this restriction can be relaxed (and indeed has been 
by Florian Bolgar in his Internship Report to the Ecole Normale Superieure 
(2012)). The background vacuum space-time with the spin vector s* = 0 
has line element given by (I96I) - (|1U2I) with F = 0 and the coefficients of 
the powers of r are given by (|104p - (jll0l) with F = 0. Thus in coordinates 
X* = {x,y,r,u), with i = 1,2, 3,4, the only coordinate component of the 
metric tensor involving the 4-acceleration a® of the world line X® = w'^{u) is 

544 = l-2r/io + 0(r2) . (127) 

The coefficient of — 2r on the right hand side here is the I = 1 spherical 
harmonic ho, which vanishes if and only if a® = 0. In this case the world line 
A® = t(;®(u) is a time-like geodesic and is thus the history of a test particle 
moving in the vacuum gravitational field modeled by the space-time with 
line element (1961) with F = 0. Now we consider a particle with spin s® 7^ 0 
having world line A® = w^{u) in the space-time with line element (j96p with 
F / 0. Now (I127p is replaced by 

544 = 1 - 2r (^ho - ^RijM F + J2) + 0 {r^) , (128) 

with J2 given by (|124l) . The first two terms in the coefficient of —2r here 
are I = 1 spherical harmonics while the third term J2 is an / = 2 spherical 
harmonic. If the equations of motion of a spinning test particle are obtained 
in the same way as those of a non-spinning test particle then equating to 
zero the I = 1 terms in the coefficient of — 2r in (I128p results in 

at F = ho = ^Rijkiiu) F s^'- , (129) 

for all possible values of F and thus the equations of motion of a spinning 
test particle, neglecting spin-spin terms, are 

ai = ^Rijkiiu)v^ , (130) 

in agreement with Mathisson [20] and Papapetrou m- 

5 Discussion 

Perhaps the most striking aspect of the technique presented here for ex¬ 
tracting equations of motion from field equations is the absence of infinities 
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arising in the calculations, either in the form of infinite self energy of the 
Dirac type [18] or in the form of divergent integrals. We have already indi¬ 
cated following ()84l) above that in the case of the DeWitt and Brehme [T8| 
work this may be due to the fact that they have considered a different prob¬ 
lem to the one we have described in Section 2. Much work in recent years 
has been done under the general heading of “the self force problem” and, 
fortunately, this work is well represented in the current volume, offering the 
interested student a direct comparison with our work and the challenge of 
connecting the different points of view. 
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